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Abstract 



Calculations with a synoptic case study show that long-wave 
radiative cooling tends to reduce the available potential energy, 
especially in the upper troposphere. Synoptic-scale precipitation 
amounts resulting from destabilization of clouds by long-wave cooling 
are computed. These range up to 1.4 mm in 12 hours. This de- 
stabilizing effect may be important in explaining the nocturnal maximum 
of precipitation over the sea. It may also contribute significantly to 
cyclone development. 



1 . Introduction 



This paper is a continuation of Scientific Report No. 1 published 
earlier under this Research Project (Danard 1968). 

Calculations of Katayama (1967) show that long-wave radiation, 
absorption of solar radiation and release of latent heat are about equally 
important in the annual mean generation of available potential energy 
over the northern hemisphere. Section 2 is concerned with the first 
of these processes. Cloud-top cooling reduces temperature contrasts 
if the cloud mass is in the warm air. 

The existence of a nocturnal maximum in frequency of precipitation 
over the sea has been demonstrated by Kraus (1963). He studied the 
records of weather ships in the eastern Pacific and North Atlantic 
Oceans. With 8-10 years of data for July, he found that precipitation 
was reported in 19 per cent of all nocturnal synoptic reports compared 
to 12 per cent of the daytime reports . 

Kraus proposed that during the day, absorption of short-wave 
radiation reduces the production of liquid water in clouds, thereby 
lowering the probability of precipitation. At night, this effect is ab- 
sent. 

The nocturnal maximum in precipitation occurrence over the central 
United States in summer was studied by Bleeker and Andre (1951). 

Their explanation was that low-level differential cooling at night due 
to variations in topography would result in upward motion between the 
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Appalachians and the Rockies. 

In Section 3 it is suggested that an important factor is the de- 
stabilizing influence of long-wave cooling at night from cloud-top. The 
amount of precipitation associated with convective overturning is com- 
puted. The destabilization may also augment the large-scale vertical 
velocity. This would further increase the precipitation as well as 
accelerate cyclone development. 
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2. Effects of long-wave radiation on the available potential energy. 



Consider the thermodynamic energy equation 



at 



= - V • 7 a+ + 



RH 

C p 
P 



(1) 



where cr = -a S in Q/3p, u.' = dp/dt, and H is the rate of input of 
heat per unit mass . In (1) , let a = a + a 1 , cr = a + a 1 , etc . , where 
the bar denotes the average over an isobaric surface of area S and 
the prime the departure therefrom. Eq . (1) may then be transformed to 
(Danard 1966 or Wiin-Nielsen 1964) 
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where 



1 r P bp <* ' 2 jcj 

A= 2^ J J — dSd P 

9 P t 3 * 



(3) 



will be referred to as the available potential energy of the region (c f. 
Eq (10) of Lorenz (1955) or Eq (2.2) of Wiin-Nielsen (1964)). In (2) 
and (3) , p and p^ denote the pressures at the top and bottom of the 
region of interest. 

If a quasi-geostrophic balance exists, w' will be related to H'. 
However, attention will be restricted here simply to the final term 
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in (2), representing the generation and designated as G: 

<5 - 5 s ; ; Pb s ^'dsdp 

P P t s cr p 



( 4 ) 



In (4), H' will be identified with long-wave radiative effects com- 
puted in Section 6 of Sci. Rep. No. 1. The area S is bounded by the 
latitudes 28. 7N and 51.3N and by the longitudes 70W and 115W. Com- 
putations are performed for three layers: 900-700 mb, 700-500 mb, and 
500-300 mb. Results are given in Table 1 . 

In the lower layers, G has small, positive values. This is due in 
large part to the absence of cooling in the lower portion of the main 
cloud mass (see Figs. 8 and 9 of Sci. Rep. No. 1.) In the lowest 
layer, there is also the effect of cooling from the top of stratocumulus 
cloud in the cold air over Kansas and Oklahoma (see Fig. 9 of Sci. Rep. 
No. 1). However, in the uppermost layer, G is negative and fairly 
large in magnitude. This is due to radiation from cloud-top in the warm 
air. Thus long -wave radiation tends to reduce the temperature contrast 
in the upper troposphere, a not unexpected result. Assuming G con- 
stant and neglecting the other terms on the right side of (2) , the quo- 
tient A/ |g| represents the time required to effect a change in A 
equal in magnitude to the initial value. Whereas this time is of the 
order of weeks or longer in the lower layers, it is of the order of days 
in the highest layer. This suggests that long wave radiation from cloud- 
top may be important in the upper troposphere even for short time periods . 
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Table 1. Comparison of terms A and G . See Eqs (2) - (4). 



Term 


Units 


900-700 mb 


700-500 


500-300 


900-300 


A 


m 25 

10 ergs 


6.04 


7.36 


4.97 


18.4 


G 


in 20 -i 

10 ergs sec 


0.26 


0.08 


-1.66 


-1.32 


a/|g| 


days 


27 


104 


3.5 


16 
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3. Effects of long-wave radiation from cloud-top on precipitation rates 
The basic principles in computing the convective precipitation are 
similar to those employed by Manabe et al (1965). Consider an unstable 
atmosphere initially at rest. If the motions arising from convective 
overturning are dissipated as heat, it may be shown that 

6 J C p T dm = J6Q dm (5) 

where 6 refers to a change over a finite time interval, dm is an 
element of mass and 6 Q is the heat added per unit mass from all 
sources other than friction. In (5), the integration extends over the 
entire atmosphere. Eq (5) is derived in the Appendix. In a dry adiabatic 
transformation the sum of the potential and internal energies (represented 
by the integral on the left side of (5)) is conserved. 

Eq (5) is applied to a layer cloud extending from 1000 to 400 mb. 

The lapse rate is assumed initially pseudo-adiabatic. To simulate con- 
ditions during the night, the cloud top radiates for 12 hours. As this 
proceeds , the cloud becomes unstable . The convective overturning 
then produces precipitation and returns the lapse rate to the pseudo- 
adiabatic value. The liquid water contents of the initial and final 
states are assumed equal (i.e. , all water which condenses during this 
process falls out as precipitation). 

Applied to the above problem (5) may be written 

Po 

6T 
P4 




L ,. P ° 

dp = AQ + - J 6r dp 
5 p 4 



( 6 ) 
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where p = 400 mb, p = 1000 mb, 6T(< 0 ) and 6 r (< 0 ) are the 
4 o 

differences between final and initial temperatures and mixing ratios, 

A Q (>0) is the heat lost by radiation, and L is the latent heat of 

vaporization. The total precipitation is 

P, 



P = -f° fir & 
P 4 



(7) 



The interval between p. and p is divided into six layers each 

4 o 

100 mb thick. The integrals are then expressed as sums of contributions 
from these layers. With a known initial state, Eq (6) is thus an equa- 
tion in 12 unknowns, the final values of T and r at the mid-points 
of the six layers. Eleven more equations are needed. Five are obtained 
from the final condition that 



ae 

€ 

5p 



= 0 



where 



T P R/C 

e ' = (T+i 1 ) (-f-') p 
Cp v p y 



( 8 ) 



(9) 



is the equivalent potential temperature. The remaining six equations are 

obtained from the Clausius-Clapeyron equation. If r. , e. and T. 

ii i 

denote initial mixing ratio, vapor pressure and temperature, one obtains 

0.622 r.L 

6 r = — £ — 6T (10) 



RT. 



P - e - 



Eq (10) is applied at the six levels. 
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In computing the radiation lost from cloud-top, equations similar 
to (2) and (13) of Sci . Rep. No. 1 were derived. The level p^ was 
taken as 400 mb (cloud-top) and q^ was set equal to the saturation 
value corresponding to the temperature there. The upward flux at cloud- 
top is thus 



F 

u 




a - « u i 



an 



where F^ is the black body flux corresponding to cloud-top temperature, 
and € is the water vapor emissivity for the layer above the cloud. The 
heat lost by radiation is , per unit area , 

AQ = F u At (12) 

where At = 12 hrs . 

The method described above is applied to three cases where the 
initial state is saturated with a constant wet-bulb potential temperature 
0 w of 0, 10 and 20 C. These would represent conditions over polar, 
mid-latitude and sub-tropical oceans, respectively. Results are shown 
in Table 2 . 

The values of A Q decrease with increasing initial 0 . This 

w 

is because the decrease in (1 - € ) in Eq (11) over-compensates the 

increase in F^ . Release of latent heat tends to offset the cooling due 

to long-wave radiation. This results in fairly small values of 6 0 . 

w 

Precipitation amounts in Table 2 are to be regarded as averages over 

large areas, and not as values associated with individual convective 
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Table 2 . Net cooling (6 0 ) , heat lost by long-wave radiation ( AQ ) 

w 

and total precipitation (P) during 12-hr period for cloud extending from 
1000 to 400 mb. 



Initial 0 

w 

(deg C) 


60 

w 

(deg C) 


AQ 


P 


9 -2 

(10 ergs cm ) 


(mm) 


0 


1.2 


11.2 


0.9 


10 


00 

• 

o 


9.8 


1 .4 


20 


0.3 


6.1 


1.4 
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clouds. The amounts increase slightly with initial 0 . Note that a 

w 

cooling of only 0.3 C for = 20 C produces more precipitation than 

does a cooling of 1.2C for 0 = 0C. 

w 

Precipitation amounts given in Table 2 are fairly modest. However, 
in a quasigeostrophic numerical prediction model, the large-scale 
vertical motion varies inversely as the value assigned to the static 
stability. Thus the cloud destabilization may lead to an increase in the 
large-scale vertical velocity also. This would contribute further to the 
amount of precipitation. 
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Appendix - Derivation of Eq (5) 

Consider the equation for 3-dimensional motion 
HV -4-4-4 

_ = _ a y p _ 2 0 x V+ g + F (Al) 

where V = u i + v j + w j , 7 = i 3/3x + j 3/dy + k 3/dz and F is the 
frictional force associated with motions arising from convective over- 
turning. Multiplying Eq (Al) scalarly by V gives 

d -* -» 

^ (y- + gz) = - a v- v P + v ♦ f (A2) 



The first law of thermodynamics may be expressed as 

dT _ dQ_ da 



v dt 



dt 



dt 



(A3) 



Adding (A2) and (A3) and using the equation of continuity 



da _ r? 

- = ar -v 



(A4) 



yields the equation 



^ (^- +gz + C y T) =- a7 • P V+^ + V * P (A5) 

Now if F is any function, use of (A4) gives 

P F + V ‘ P F V (A 6) 

Integrating (A5) with respect to mass over the entire atmosphere, and 
making use of (A 6) and the relation 

00 QO 

J p(gZ + C T) dz = J* pC T dz (A 7) 

o o p 

one obtains 



11 



afl + V )dm = ^ dm + J ' 7 ' ?dm <*<» 

where dm = pdz dA is an element of mass. 

Now if the motions arising from convective overturning are dissipated 
as heat, 

v 2 

6 J — dm = 0 (A9) 

Furthermore , 

J 5 Q^dm = -J[J > V*Fdm]dt (A 10) 

where 6 is the heat added per unit mass due to friction and the time 
integration extends over the interval to which refers . Thus (A8) 

becomes 

6 J CpT dm = J 6 Q dm (5) 

where 6 Q refers to heat added from sources other than friction. 
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